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c3 . Abstract 

' Magnetic excitations in an array of (VO)2P207 single crystals have been 
g' 

I ' measured using inelastic neutron scattering. Until now, (VO)2P207 has been 

^ \ thought of as a two-leg antiferromagnetic Heisenberg spin ladder with chains 

^ ' running in the a-direction. The present results show unequivocally that 

^ _ (VO)2P207 is best described as an alternating spin-chain directed along the 

00 

' crystallographic 6-direction. In addition to the expected magnon with mag- 

o ■ 

[ netic zone-center energy gap A = 3.1 meV, a second excitation is observed 

t — 

I at an energy just below 2A. The higher mode may be a triplet two-magnon 

' bound state. Numerical results in support of bound modes are presented. 
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The S=l/2 alternating Heisenberg chain (AHC) is a fascinating quantum system that is 
currently the subject of much interest. We have established that the material (VO)2P207 
, previously considered to be a spin ladder, is in fact an excellent realization of the AHC. The 
physics of the AHC is also very relevant to spin-Peierls materials such as CuGeOS 0. Recent 
theoretical work on the AHC underscored the potential importance of two-magnon bound 
modes. In this paper, we review our neutron scattering experiments on the alternating chain 
material (VO)2P207 . In addition to the expected one-magnon excitations, we observe an 
extra mode, which may be a two-magnon bound state. We follow with a discussion of some 
related theoretical issues. 

The crystal structure of (VO)2P207 is nearly orthorhombic, with a slight monoclinic 
distortion so that the space group is P2i 0. The room temperature lattice parameters 
are a=7.73A, 6=16.59A, c=9.58Aand /5=89.98°. The magnetic properties of (VO)2P207 
arise from S = 1/2 V^^ ions situated within distorted VOe octahedra. Face-sharing pairs 
of VOg octahedra are stacked in two-leg ladder structures oriented along the a-axis. The 
ladders are separated by large, covalently bonded PO4 complexes. The structure is illustrated 
schematically in figure 1. 

The susceptibility of (VO)2P207 powder can be accurately reproduced by either a 
spin ladder (with J\\ ^ J±) or by an alternating chain |^^, but the expectation that the 
PO4 group would provide a weak superexchange path led to a general acceptance of the 
spin ladder interpretation of (VO)2P207 . Pulsed inelastic neutron scattering measurements 
on (VO)2P207 powders showed a spin gap of 3.7 meV, which was interpreted as further 
support for the ladder model. 

Beltran-Porter et al. examined the superexchange pathways in several vanadyl phos- 
phate compounds, and were led to question the spin-ladder interpretation of (VO)2P207 
. Instead they proposed that an alternating V-O-V-PO4-V chain in the 6-direction was a 
more likely magnetic system. The observation of a second spin excitation near 6 meV (not 
predicted by the ladder model) in a recent triple-axis neutron scattering experiment on 
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(VO)2P207 powder, and the discovery of strong superexchange through PO4 groups in the 
precursor compound VODPO4 ■ I/2D2O ||T0|, also cast doubt on the spin ladder interpre- 
tation. For these reasons we undertook studies of the spin dynamics in (VO)2P207 single 
crystals. 

Measurements of the excitations were made using a single crystal array of approx. 200 
oriented (VO)2P207 crystals of typical size 1x1x0.25 mm^. The resulting sample had an 
effective mosaic spread of 8 — 10° FWHM. Inelastic neutron scattering measurements were 
carried out using triple-axis spectrometers at the HFIR reactor, Oak Ridge National Labo- 
ratory; full experimental details can be found in |l[]. 

Scans [HP at T = lOK showed two modes of roughly equal strength at the antiferromag- 
netic zone-centre (0, vr, 0) at energies of = 3.12(3) meV and A„ = 5.75(2) meV. Full 
resolution convolutions with the fitted dispersion showed these modes to be resolution lim- 
ited. The disappearance of both modes at higher temperatures confirmed their magnetic 
origin. The modes were found to track approximately in Q close to (0, vr, 0) (see Figure (2)). 
At the zone-boundary (0, 7r/2, 0) only a single mode was observed at an energy of ~ 15meV. 
Because of poor instrumental resolution it was not possible to tell whether the modes had 
coalesced or were simply not resolved. 

Fig. 2 shows the measured dispersion for both modes along a*, h* and c*. The excitation 
energy is almost independent of Qc (middle panel), implying a very weak coupling along 
c. The dependence of energy on Qa is much weaker than on and is ferromagnetic. The 
strong Qfe dependence implies that the exchange coupling is dominantly one-dimensional 
along the 6-direction, confirming the V-O-V-PO4-V alternating chain proposed in |^ and 

For any exchange alternation — as occurs with the two inequivalent exchanges along h — 
a gap should appear in the dispersion (as observed) and the absence of magnetic ordering in 
(VO)2P207 is consistent with a singlet ground state. However the observation of an extra 
mode requires a more thorough theoretical investigation. Since the high temperature limit 
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of the magnetic susceptibility is consistent with expectations for simple S = 1/2, g = 2 
spins, the possibility that the upper mode is an additional low lying single ion excitation 
can be ruled out. Two other plausible explanations of the second peak are (a) splitting due 
to an exchange anisotropy, and (b) a physical two-magnon bound state. 

Although a pseudo-Boson calculation including exchange anisotropy gave an excellent fit 
to the dispersion — the solid line in Figure (2) is a fit to this model [|l| — considerable 
exchange anisotropy (~ 15%) was necessary to account for the mode splitting. Recent 
single crystal magnetic susceptibility measurements [Til were quantitatively consistent with 



the previous powder results |^ and found little if any evidence for anisotropy. Also the 
coupling in the precursor compound VODPO4 • I/2D2O was found to be consistent with 



isotropic exchange |jTO| suggesting that one should seek another explanation for the second 
mode, and because the energy of the upper mode at (0,7r,0), A^j, is just below 2/S.i, a 
bound two-magnon mode may provide a good explanation. In support of this explanation, 
preliminary high-field measurements show splitting of both modes which is consistent with 
both modes being triplets []T2[. 

To gain some insight into the formation of two-magnon modes in (VO)2P207 we have 
studied the = 1/2 AHC using numerical techniques. The isotropic AHC Hamiltonian is 

L/2 

H = ^ J 5*24-1 • + aJ 821 ■ 821+1 , (1) 
1=1 

where J > and 1 > a > 0. Equation (1) has been studied analytically and numerically 
over many years, but it had not been appreciated until recently that 8 = and 5 = 1 bound 
magnon states may form. Uhrig and Schulz have used field theory and RPA methods to 
study these modes at = 7r/2 and k = 0,71. The existence of these bound states depends 
subtly on the kinetic and potential energies of pair formation, and occur for only certain 
values of k. 

Perturbation theory in a about the dimer limit {a = 0) provides a quantitative basis for 
understanding the excitations for small a, and also provides insight into the competition 
between potential and kinetic energy effects in bound states [|1^]. Figure (3) shows the one- 



and two-magnon excitation spectra calculated within a simplified approximate first order 
(one- and two-magnon manifold) treatment of the AHC At k = 7r/2 there is a node in 
the two-magnon continuum which corresponds to a degeneracy in the total kinetic energy 
uj{ki) + 00(71/2 — ki) of two magnons. The S = and S = 1 bound states lie well below the 
continuum lower boundary. However at = and vr only the S = bound state is seen. The 
continuum is much broader at A; = 0, vr indicating larger mixing effects which disrupt the 
S = 1 bound state. Although no S* = 1 bound state forms, the attractive potential still leads 
to a strongly enhanced scattering cross-section S{Q,uj) at the continuum lower boundary 
1131 , see dashed line in Figure (4). The 5=1 bound state appears clearly at the k = ti/2 
point (solid line in Figure (4)). It should be noted that the neutron scattering cross-section 



for the S = Q mode is zero, however this mode may be visible by light scattering [13 



Harris [|T4| used a reciprocal space perturbation theory to calculate the ground state and 
excited state energy up to Oi^or'). This gives a = 0, vr energy gap of Egap = J(l — a/2 — 
3q;^/8 -|- a^/32). However these results can be derived more easily in real space [Q, and in 
the case of the ground state energy, we have extended the calculation to O(a^), 

eo{a)/J = -3/2^ - (3/2^) • - (3/2«) • 
-(13/2^2) . ^4 _ (g5/3) . (^/2i4) ■ a' - 0{a''). (2) 

The perturbation series appears to be rapidly converging for a < 0.5, and may have a radius 
of convergence of unity. 

Because (VO)2P207 has a ~ 0.8 [|T|, we have used a numerical Lanczos algorithm on 
finite L = 4n lattices of up to L = 28 and with approximately 14 place accuracy to study 
the ground states and binding energies up to similar values of a. Full details will be given 



elsewhere [0]. Figure (4) shows the calculated binding energies of the S = bound mode 
at k = tt/2, and A; = 0, as well as those for the S = 1 bound mode at k = tt/2. The results 
show strong binding at a = 0.8 of the S = mode at A; = 7r/2 but the situation is not clear 
for A; = 0. They also suggest weak binding for the S=l mode at tt/2 at the alternation for 
(VO)2P207 . Unfortunately finite size effects precluded an accurate determination of this 



binding energy. 

In order to make a quantitative comparison with (VO)2P207 S{Q,uj) is required for the 
bound modes and continuum. We are currently undertaking calculations to quantify this. 
The effects of interchain couphng have been neglected and these may enhance the binding. 
Next-nearest neighour exchange within the chains may have a similar effect. The a pertur- 
bation theory provides a useful quantitative guide to such effects, and further theoretical 
studies are in progress. We also note that similar dynamics are also important in many 
other low-dimensional Hamiltonians such as spin ladders and we shall present some work on 
those in the future. 

In conclusion, we have measured an extra mode in the alternating chain system 
(VO)2P207 . The evidence suggests that this is a two-magnon bound state. Perturbation 
theory and Lanczos calculations give an insight into the formation of bound modes. 

We thank J.Thompson for sharing his susceptibility results with us prior to publication. 
Oak Ridge National Laboratory is managed for the U.S. D.O.E. by Lockheed Martin Energy 
Research Corporation under contract DE-AC05-96OR22464. Work at U.F. is supported by 
the U.S. D.O.E. under contract DE-FG05-96ER45280. 
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Figure Captions 



Figure 1. Schematic depiction of the structure and magnetic interactions in VOPO. 
The spin ladder model previously thought to describe VOPO has nearest neighbor exchange 
constants J\\ along the a ( "ladder" ) direction and J± along the h ( "rung" ) direction. In the 
alternating chain model, nearest neighbor ions are alternately coupled by constants Ji 
and J2 along the h (chain) direction. Neighboring spins in adjacent chains are coupled by 
J a- Magnetic coupling in the c direction is negligible. 

Figure 2. Measured dispersion of magnetic excitations in VOPO at T = lOK. When not 
visible error bars are smaller than the size of the plotted symbols. Filled circles (open dia- 
monds) are points from the lower (upper) energy mode. The solid lines are dispersion curves 
calculated using parameters obtained by fitting to a pseudo-Boson model [^. Wavevectors 
are plotted in units corresponding to the VOPO reciprocal lattice. 

Figure 3. Schematic depiction of the one- and two-magnon excitation spectra of the 
S = 1/2 AHC with an alternation of a = 0.2. An S = 1 bound mode appears below the 
continuum at ~ 7r/2. The more deeply bound 5 = mode (dashed line) is not visible to 
neutrons scattering. 

Figure 4. Calculated S{Q^uj) for constant-Q scans aX k = 7r/2 (solid line) and k = ir 
(dashed line) using the first order perturbation approach with a = 0.2. 

Figure 5. Calculated binding energies of the S = and S = 1 bound states using a 
Lanczos method |TB| at A; = and 7r/2. The binding energies are given in units of J. 
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